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Abstract 


The wide variety of available interactive methods brings the need for cre- 
ating general interactive algorithms enabling the decision maker (DM) to 
apply freely several convenient methods which best fit his/her preferences. 
To this end, in this paper, we propose a general scalarizing problem for multi- 
objective programming problems. The relation between optimal solutions of 
the introduced scalarizing problem and (weakly) efficient as well as properly 
efficient solutions of the main multiobjective optimization problem (MOP) is 
discussed. It is shown that some of the scalarizing problems used in different 
interactive methods can be obtained from proposed formulation by selecting 
suitable transformations. Based on the suggested scalarizing problem, we 
propose a general interactive algorithm (GIA) that enables the DM to spec- 
ify his/her preferences in six different ways with capability to change his/her 
preferences any time during the iterations of the algorithm. Finally, a numer- 
ical example demonstrating the applicability of the algorithm is provided. 
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1 Introduction 


The general goal of solving a multiobjective optimization problem (MOP) is 
to support the decision maker (DM) seeking the most preferred solution of 
many Pareto optimal solutions as the final one. Inasmuch as finding a most 
preferred solution needs some extra information from the DM, interactive 
approaches, based on the participation of the DM, have become popular. 

In interactive methods, an iterative algorithm is proposed. Then, the 
steps of the algorithm are repeated where at each iteration, some informa- 
tion is given to the DM and he/she specifies his/her preferences. The process 
is repeated until the DM is satisfied with regard to the obtained solution. 
The benefits of using interactive approaches are that, the DM (2) does not 
need to have any global preference structure, (ii) has the possibility of learn- 
ing about the interrelationship between the objectives, (<i) can learn about 
the feasibility of solutions during the solving process. 

Heretofore, many interactive methods have been suggested in the litera- 
ture [1,13, 19, 23, 26, 30,31]. As pointed out already, interactive methods are 
very useful and realistic to solve an MOP. However, since there have been 
many interactive methods available, it is not easy to choose an appropriate 
method conveniently. Therefore, creating global algorithms with an ability 
to accommodate different methods will be useful. By creating a global algo- 
rithm, it is possible for the DM to select freely an appropriate method (and 
the way of specifying preference information) as well as to switch between 
methods. To this end, it is necessary to design a general scalarizing problem 
yielding scalarizing problems used in different interactive methods. 

Until now, some global algorithms have been proposed. For example, 
Gardiner and Steuer [7,8] proposed a unified algorithm including nine to 
thirteen different methods. Romero [27] presented another general optimiza- 
tion structure, called extended lexicographic goal programming. Moreover, 
Vassileva [32] suggested a general scalarizing problem which incorporates 
different scalarizing problems. More recently, based on a global formulation 
(GLIDE), Luque et al. [21] proposed a global procedure which accommodates 
eight interactive methods of different types. Nevertheless, their formulation 
is unlikely to consider the computational efficiency, therefore Ruiz et al. [28] 
improved the computational efficiency of GLIDE by reformulating it. 

In some of the mentioned publications, the authors have provided theo- 
rems concerning (weak) efficiency of the optimal solutions of their proposed 
general scalarizing problems [21, 28,32] and as far as we know few results 
related to proper efficiency have been provided. Now, in this paper we sug- 
gest a general scalarizing problem which not only considers computational 
efficiency by reducing the number of added constraints, but also provides 
theorems concerning (weak) efficiency as well as proper efficiency of its op- 
timal solutions. The provided results are established without any convexity 
assumption. Also, by setting suitable values for parameters and index sets of 
the proposed general scalarizing problem, we obtain many known scalarizing 
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problems. Based on the mentioned problem, we propose a general interactive 
algorithm (GIA) to solve a given MOP, subsequently. In this algorithm, the 
DM has the ability to specify his/her preference information in six different 
ways. 

The rest of this paper is organized as follows: Section 2 contains some 
preliminaries and basic definitions. In Section 3, we propose our general for- 
mulation and obtain some theorems. Section 4 gives some scalarizing prob- 
lems used in different interactive methods which can be obtained from our 
general formulation. Section 5 contains our proposed interactive algorithm. 
In Section 6, some computational and theoretical advantages are mentioned. 
An example is presented in Section 7 and finally, in Section 8 conclusions are 
given. 


2 Preliminaries and basic definitions 


A general multiobjective optimization problem can be written as: 


(MOP) min f(x) (1) 
st. XEN, 
where ¥ C R” isa nonempty compact set, and f(x) = (fi(x), fo(x),..-, fp(x)) : 
X — R? is a vector-valued function. 

The set of all attainable outcomes or objective vectors is defined as the 
image of the feasible solutions x € ¥ under f. In fact Y := f(¥) C R?. 
For y? and y? € R?, y* < y? means that yj < y?, for each i = 1,--- ,p, 
also yt < y? stands for yt < y? and yt ¥ y?. Furthermore, y! < y? 
means that y} < y?, for each i = 1,...,p. The Pareto cone is defined as 
RS ={fyeR’:y2 0}. RE and R& are defined, similarly. In this paper, we 
shall assume that Y := f(4) is bounded. 

Definition 1. A feasible solution & € X is called: 


(i) weakly efficient (weakly Pareto optimal) solution to MOP (1) if there 
is no other x € X such that f(x) < f(&), 


(it) efficient (Pareto optimal) solution to MOP (1) if there is no other x € 
X such that f(x) < f(&), 


(itt) properly efficient (properly Pareto optimal) solution to MOP (1) if it is 
efficient and there exists a real positive number M such that for each 
i € {1,2,...,p} and each x € & satisfying fi(x) < fi(x), there exists 
an index j € {1,2,...,p} with f;(X) < f;(x) and 


Fi(X) — fil) 


file) = f)(@) = 
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The set of all weakly efficient, efficient, and properly efficient solutions 
of MOP (1) will be denoted by we, Vz and pz respectively. The image 
f(x) € Y of an (weakly, properly) efficient solution x € 4, is called (weakly, 
properly) nondominated point. 


Remark 1. Obviously, ¥pz C Vz C Xwe. 


Remark 2. In this paper, we use definition of proper efficiency in the sense 
of Geoffrion [9]. There are other definitions of proper efficiency which are 
almost the same when using the Pareto cone as the order cone. For consid- 
ering relationships between different definitions of proper efficiency one can 
refer to [4]. 


Definition 2. The ideal point y! = (yf,... 02) of MOP (1) is defined by 
ys = Mingex fi(x), i = 1, ott yp. 


Definition 3. The point yU := y!— a, where a € R& is a vector with small 
positive components, is called the utopia point of MOP (1). 


Definition 4. The nadir point yN = (yN,... yp) of MOP (1) is defined by 


ye = MALxE Xz fi(X), 1 = 1, wth yp. 


Definition 5. The vector ¥ = (¥1,-.-.,Y%p) € R”, consisting of the desired or 
aspiration values to the DM, is called a reference point. It should be noted 
that reference point may be achievable or not. 


One of the most popular approaches to solve a given MOP is scalariza- 
tion, which involves formulating a single objective problem associated with 
the given MOP. Let us consider a single objective programming problem as 
follows: 


min g(x) 
s.t. x ES, 


(2) 


where g: SR. 


Definition 6. A feasible solution x € S is said to be 
(i) an optimal solution of problem (2) if g(&) < g(x) for allx ES, 
(ii) a strictly optimal solution of problem (2) if g(X) < g(x) for allx € S\{x}. 


3 A general scalarizing problem 


In this section, we propose a general scalarizing problem associated with MOP 
(1), which is defined such that many scalarizing problems, used in different 
interactive methods, can be deduced from it by selecting suitable values of 
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parameters and index sets. The general scalarizing problem is proposed as 
follows: 


P 
i Bt Fy are : x)— rk 
min pas Xj (fil ) ee v0 Ns we (fe(X) ) 


st. 4 oe 


(3) 


where \F > 0, p > 0, 6%, r¥, and wf > 0 are parameters specified depending 
on the information given by DM. Also, I 4 0 and I* are index sets, which 
are subsets of {1,--- ,p}. Notice that hereafter, we make the assumption that 
on the proposed scalarizing problem, the parameters 5*,i € I* are selected 
such that problem (3) remains feasible. Let & be the current iteration. Then, 
the optimal solution obtained from scalarizing problem (3) is defined by x*++ 
and the corresponding objective vector by f(x*+!). 


According to [14, p. 305], or [23, p. 97], if we replace the max term by a 
new variable z € R, then problem (3) is equivalent to the following scalarizing 
optimization problem: 


min z 


Pp 
AF (fal) —rh tp wh (fix) —rk)) <z vielr, 


t=1 
filx) < oF Vi € IF, 
xENX. 


S.t. 


Notice that the scalarizing problem (3) is nondifferentiable, even if the main 
MOP (1) is differentiable (i.e., the objective functions and constraint func- 
tions are differentiable). Therefore, if the original MOP is the differentiable 
we propose to use formulation (4) since it preserves differentiability. In this 
case, the scalar optimization problem (4) can be solved with standard meth- 
ods of (non)linear constraint optimization or using available single objective 
solvers. However, if the original MOP (1) is nondifferentiable, both scalarized 
problems (3) and (4) are nondifferentiable, too. In this case, the scalarized 
problem (3) is recommended since it has a reduced number of constraints. 


It should be noted that, unlike the formulations proposed in [21,28], the 
bounds on trade-offs generated by the suggested formulation are independent 
of parameters A;. For more details about bounds on trade-offs see [17, 18]. 
So far, many authors have provided theorems concerning weak efficiency and 
efficiency of the optimal solutions of the scalarized problems used in the 
interactive methods. Now, we prove some general theorems concerning weak 
efficiency, efficiency, as well as proper efficiency of (strictly) optimal solutions 
of problems (3) and (4). It is important to point out that the following 
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theorems are general and many theorems concerning (weak, proper) efficiency 
[4,23] can be resulted from them. Moreover, the theorems are provided with 
no convexity assumption. Since problems (3) and (4) are equivalent, we only 
provide theorems for the first one. 


Theorem 3. Let \* > 0 Vi € Ik. If +! © X is an optimal solution of 
problem (3), then x**+" is a weakly efficient solution of MOP (1). 


Proof. Let x*++ € X be an optimal solution of problem (3) and suppose that 
x*+1 ¢ Xyw. Then, there exists x € XY such that f(x) < f(x**t+). Therefore, 
filx) < fi(e*t*) < oF Vi € IX, which means x € & is a feasible solution for 
problem (3). Also, we have 


and 
Pp Pp 
p>. wt (felx) — rf) < p >> wh (fe(R**") — rf). 
t=1 t=1 
Therefore, 
k k41 k ek k+1 k 
okt okt 
Des (file yar; + ed we hele )-7h)) > 
Pp 
are (fil) —re y+ pd, we hela) > A), 


which is a contradiction with optimality of @*+!. Thus, x*+! € Xwe. 


In the following theorem, utilizing the general formulation (3), a sufficient 
condition for efficiency is provided. 


Theorem 4. If x**+! € X is a strictly optimal solution of problem (3), then 
xh € XR. 


Proof. The proof is similar to that of Theorem 3. 


It is found out from part (iz) of Definition 1, that in an efficient solu- 
tion it is not possible to improve any criterion without deterioration of at 
least one other criterion. Sometimes, these trade-offs may be unbounded and 
it is obvious that efficient solutions with bounded trade-offs (called prop- 
erly efficient) are desirable. Until now, many scholars have considered rela- 
tionships between optimal solutions of the scalarizing problem used in their 
proposed interactive methods and (weakly) efficient solutions of the related 
MOP [21, 28,32], but there are fewer results concerning proper efficiency. In 
the following theorem, we provide a sufficient condition concerning properly 
efficient solutions of MOP (1). 
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Theorem 5. If x*+! € X is an optimal solution for problem (3) with * > 
OVie TT, p>O and wk ERS, then x**1 € Xpp 


Proof. We show that #*+! € Vp. Let x*+! ¢ Xp. Then, there exists x € X 
with fi(x) < fi(&*t"), Vi © {1,--- ,p} and f;(x) < f;(x**1) for some j € 
{1,--- ,p}. Hence, fi(x) < fi(x i) < 6% Vi € IX. Thus, x € & is a feasible 
solution of (3). Using the assumptions and the definition of efficiency, it 
follows that: 


max \¥ (file Rly, age i p>o uk (fe( (a*+1) — rf)) > 


ie rk t=1 


max Af ( (x) —r? + w —T i 
sere Fil) edu i (f(x ‘) 
This is a contradiction with optimality of #*+1 and therefore #*+! € Xp. 
Now, we show that #**1 is a properly efficient solution to MOP (1). To this 
end, we define: 


and consider an index i € {1,--- ,p} and X € ¥ such that f;(X) < f;(#**?). 
To prove the proper send of <*+1, we must show that there exists an 
index j € {1,2,...,p} with f;(x**1) < f;(X) such that 


fiR"**) — fil) 


= : <M. 
Fj(%) — f;(®**") 
From efficiency of #**!, we conclude that there exists an index t € {1,--- ,p} 
such that f,(x**1) < ful ). We define 
Fi(R°) — F(X) = min (fi (&**") — fin (X)). (5) 
me{l,-- ,p} 


It is obvious that f;(x*t1) — f;(X) < 0. 
Moreover, optimality of ¢**+1 for problem (3), concludes 


max Am (fin (& jek, +0 atl filX )-r)) > 


merr 

Pp 
max Matinee) - re + p>_ wh (fi(x*™) = 7) ye 
melt t=1 


Now, let 


N (fil) 7 K+ p Souk filX )-rh)) = 


t=1 
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pr m( fn rare +o Da (fe(& )-rh)). 


Hence, 
(fil) a + pS out ess )-rh)) 2 
t=1 
P 
mars Am (Jma(°*1) — rn tO wh A) — rf) > 
P 
Ni (FiR**) — rf + pO wk (fle) — rf). 
t=1 
Then, 
0> (f(x***) m)) +0 wth xT) — f(X)). (6) 
Now, from (5) and (6), we have: 
0= (fi (% a) x) +0 wth sai — fi(X)). 


That is, 


pk (fiR**) — ful) < FR) — FR) + p Do wk (f®) — AR) < 
Zi 


(L+ p>) wh)(fy(%) — fi(@**)). 


t=1 
thi 
Hence 7 ‘ 
1 - 
F(R") - fil) - pditat we 
: < <M, 
fi) — f(T) ~ pwr i 


which completes the proof. 


It should be noted that, using suitable values for parameters in (3), we can 
provide necessary conditions related to (weakly, properly) efficient solutions 
of MOP (1). For example, if we choose IF = {1,--- ,p}, IF = 0, r®¥ = y¥ 
and wk =1 Vi € {1,...,p}, then we have the modified weighted Tchebycheff 
method [15] and, by Theorem 4.2 in [16], for every properly efficient solution 
of MOP (1) we can find suitable parameters \¥ > 0, Vi € {1,---,p} and 
p > 0 such that this properly efficient solution be an optimal solution of (3). 
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4 Achieving different scalarizing problems from the 
general formulation 


The general formulations (3) and (4) are generalizations of already known 
scalarizing problems. In this section, we are going to show that how many 
famous scalarizing problems (used in different interactive methods) can be 
attained from (3) and (4) by choosing appropriate values of parameters and 
index sets. We obtain the scalarizing problems from (3). By a similar method 
it is possible to obtain them from (4). 


4.1 GUESS method and STOM 


The GUESS method is one of the simplest interactive methods, proposed by 
Buchanan [2]. In this method, the DM has to determine the components of 
the reference point (y!) as preference information. At the kth iteration, the 
scalarizing problem used in this method is formulated as follows: 


(xy) — gk 
min max Aix) — i 
ae, | Yi — 9; (7) 
st. xEX 


Notice that the reference vector specified by the DM, must be strictly lower 
than the nadir objective vector, that is, y < yN. This scalarizing problem 
can be achieved from (3) by considering the following replacements: 


(1) Ik = {1,...,p} and If = 0; 


(2) wt =0, 4¥ = phe, and p= 0; 


(3) 7% = of and t= 1,205: 


The satisficing trade-off method (STOM) [24] can be obtained from (3), simi- 
lar to the GUESS method, by setting A¥ = ae and r? = yY (i =1,...,p). 
Other parameter values and index sets are the same as those of GUESS 
method. In this method, ¥ must be chosen such that y > y¥. 


4.2 Reference direction approach 


In this method, a vector from the current iteration point to the reference 
point (a reference direction) is projected onto the efficient set [20]. To obtain 
the points along the reference direction at the kth iteration, the following 
scalarizing problem needs to be solved: 
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fi(x) — (ff + td?) 
i=1,...,p Ly (8) 
St. xEX, 


where, f* is the current nondominated objective vector, d* = ¥* — f*, t has 
different discrete nonnegative values, and yp is a weighting vector that can be 
either a reference point presented by the DM or defined as yN — y¥. This 
problem can be obtained from (3) by considering the following replacements: 


(1) Ik = {1,...,p} and If = 0; 


(2) wh =0, AF = jp and p= 0; 


(3) ré¥ = fk +td* andi =1,...,p. 


4.3 Step method 


The step method is one of the first known interactive methods [1]. Eschenauer 
et al. [6] extended this method to nonlinear problems. In this method, based 
on the current objective vector (f*), the DM can improve some unacceptable 
objective functions f; (i € Jt) by relaxing some other objective function(s) 
fi @ € J&) such that JP U JF = {1,...,p}. In this regard, the DM must 
specify upper bounds «¥ > f* for functions f; (i € J#). In this case, the 


scalarizing problem is formulated as follows: 


j=l "I 
u(x) < fF Vie JE, (9) 
s.t. fix) <e* Vie Jk, 
XEN, 
I 
where e; = 4 (He ), i=1,...,p (the denominators are not allowed to be 


zero). We can obtain (9) from (3) using the following replacements: 
(1) Fh 1h cc ands es 
(2) wk =0, Vie {1,...,p}, AF = Se Vi € {1,...,p} and p = 0; 
(3) r§ =yf, Vie {1,...,p}; 


(4) 5% = fF fori ce JF and 6% =e# forie J. 
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4.4 SPOT method 


In the SPOT method, given the current objective vector f*, the DM is asked 
to select a reference objective function f; and then compare each objective 
function f; (¢ = 1,...,p, « # 1) with f; by providing the marginal rates 
of substitutions (MRSs) mi (i = 1,...,p, i # 1) [29]. The MRSs can be 


7 
approximated as mk ~ ate, i = 1,...,p, where Af? is the amount of 


improvement, provided by the DM, on the value of the objective function f; 
that can exactly compensate for the given amount Aff to be deteriorated 
of the reference objective f;. The intermediate single objective optimization 
problem, used in this method, can be formulated as follows: 


min f(x) 


“ )xen, 


where pi.,i # 1 are K.K.T multipliers, corresponding to the current non- 
dominated objective vector [21] and several values for a are set and in this 
way different solutions are obtained. This problem is achieved from (3), by 
considering the following transformations: 


(1) If = {0} and If = {1,...,p}\{4}; 
Oy Sp = ond 0) 


(3) OF = fF +a(uk —mi),i=1,...,pandiAl. 


4.5 Modified reference point method 


This method is an interactive reference direction method for solving convex 
nonlinear integer problems [31]. Here, the DM is asked to set his/her pref- 
erences as aspiration levels of the objective functions at each iteration. Let 
Jt be the set of indices of the objective functions which the DM wants to 
improve and J} denotes the set of indices which can worsen and de contains 
the indices that are satisfactory to the DM. The scalarizing problem used in 
this method is formulated as follows: 


fix) - oF Gd = Hy 
fg" eae 
- oe Vie Jy, 


xen, 


min max { 
ieJE, jet 
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where the denominators must be positive. By the following replacements, 
problem (11) can be resulted from (3): 

(Ly = Pe and 1) =F 
(2) ire A = Gre g Vie Jf, f= 
=1,...,p; 
(3) rF =y' Vie JR, rF = fF Vie JR and 6% = fF Vie JF. 


gape Vi € J}, p = 0 and 


4.6 RD method 


The reference direction (RD) method was proposed in [25]. At the kth itera- 
tion, the DM is asked to specify a reference point y”. Specifying a reference 
point is equivalent to classifying the objective functions in three classes Jf , 
J& and J*, where these index sets are the same as those defined before. The 
scalarizing problem related to the RD method is as follows: 


k 
minmas 
filx) < ft Vie JE, (12) 
st. § fi(x) <P +a(fh—g) Wie JY, 
xen, 


where 0 < @ < 1 and the denominators must be positive. The general 
formulation (3) can be transformed to RD problem (12) by the following 
replacements: 

(Ly Th SP anid Th ee 


(2) aera = Vie JF and p=0; 


Fee uF 
3 = fivie JE, oF = fF Vie JE and 6% = yg? + a( fF — gF) Vi € JE. 
3 


4.7 e—Constraint method 


In this method, one of the objective functions is minimized, while the other 
objectives are transformed into constraints by setting an upper bound [4,23]. 
The problem to be solved has the following form: 


min f;(x) 
es fo Vre tle £h (13) 


xEX. 
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By the following transformations, problem (13) can be attained from (3): 


(1) Tf = {0} and I} = {1,...,p}\{0}; 
(QA ST, p= OL en, oF ek, j 1,...,pand jl. 


4.8 The weighted sum method 


In this method, a weighting coefficient is associated with each objective func- 
tion and then the weighted sum of the objectives is minimized [4,23]. Ac- 
cordingly, solutions are obtained by solving the following problem: 


P 
p k 
min Soph fi(x) 
~ an 
st. xe XX, 


with w® > OVi € {1,...,p} and 7?_, wp? = 1. By the following replacements, 
we can obtain this problem from (3): 


(1) If = {I}, where 1 is an index with pi #0 and If = 0; 


(2) = 4, p= 4, uF pe VAL, wh HE and rk = 0 Vi € {1,...,p}. 


4.9 Hybrid method 


The hybrid method is a combination of the weighted sum method and the 
€—constraint method [4,23]. This problem has the following form: 


P 
min bf) 
aa F (15) 
s.t ix) Se, ¥j € {1,...,p}, 
| x eX, 
where pf > 0 Vi, S7P_, uP = 1 and e* = (ef,...,e%) is an upper bound 


vector. One can find this problem from (3) by the following transformations: 
(1) I* = {1}, where J is an index with yp? #0 and I¥ = {1,..., p}; 


k k 
(2) MM =D, p= Ge, wh = wy Vi A Land wf = 4; 


(3) r® =0 and 6F =e*, Vi € {1,...,p}. 


92 M. Ghaznavi, M. Ilati and E. Khorram 


Remark 3. Using similar procedure, we can obtain some other single ob- 
jective problems used in different interactive approaches. For example, the 
intermediate problems of the interactive surrogate worth trade-off (ISWT) 
method [3] and the PROJECT method [22] can be obtained easily from our 
formulations. In addition, the weighted Tchebycheff scalarizing problem [4] 
and the modified weighted Tchebycheff problem [15] are resulted from the 
proposed general scalarizing problem. 


5 General interactive algorithm 


Based on the general formulations given in Section 3, we present a general 
interactive algorithm (GIA). The proposed algorithm allows the DM to spec- 
ify his/her preference information in six different ways. Moreover, he/she 
will be able to change his/her preference information in each iteration. In 
addition to widely used ways (reference point specification, classification of 
the objective functions, and specification of marginal rate of substitution) for 
specifying preference information ( [21,28]), GIA allows the DM to specify 
his/her preferences as criteria weights, ¢«—constraint (choosing a reference ob- 
jective function and setting upper bounds for the other objective functions), 
or criteria weights and upper bounds for objective functions, simultaneously. 
The main steps of the GIA are given in Algorithm 1. 

As pointed out in Step 4, the values of parameters and index sets depend 
on the type of preference information given by DM in Step 3. For example, if 
DM specifies his/her preferences as reference point, we should set parameters 
and index sets in (3) or (4) so that one of the reference based on scalarizing 
points problems (see, for example, (7) and (8)) be attained. 


6 Computational and theoretical advantages 


The GIA and the proposed scalarizing formulation has a number of potential 
advantages both in theoretical and computational points of view. Here, we 
indicate only some key potential advantages, with special attention to those 
not shared by other competing algorithms. 


(a) Taking the special characteristics of the problem into account, GIA 
allows using more efficient optimization methods. In the first step of 
the GIA, the type of problem (differentiable or nondifferentiable) is 
specified. This step provides some advantages. For example, if the 
proposed problem is differentiable, corresponding scalarizing formula- 
tion preserves differentiability and can be solved using available single 
objective solvers. On the other hand, for a nondifferentiable problem, 


An interactive algorithm for solving multiobjective optimization problems ... 


Algorithm 1. General interactive algorithm (GIA) 


Step 1- Determine type of the MOP being solved (differentiable or nondifferentiable). 
Step 2- Compute ideal and nadir points. Set k = 1. Determine an initial solution (can be 


specified by the DM or by solving an arbitrary scalarizing problem). Denote this 
initial solution by x* and corresponding objective vector by f(x*). If the DM is 
satisfied with this solution, go to Step 6. 


Step 3- Ask the DM to provide his/her preference information based on f(%*). The DM 


can specify his/her preference information in one of the following ways: 


3.1. Specifying the desired objective function values as components of the reference 
point (gf, i=1,...,p); 


3.2. Classifying the objective functions into two classes Je and iis or three classes 
Jt ; Je and JE, described in the text; 


3.3. Specifying the marginal rates of substitutions (MRSs); 
3.4. Determining the criteria weights; 
3.5. Providing preferences with the help of e—constraint; 


3.6. Defining preferences with the help of criteria weights and selecting the upper 
bounds for all objective functions, simultaneously. 


Step 4- Based on the preference information, given by the DM in Step 3, set appropriate 


values for parameters and index sets in formulation (3) (for nondifferentiable MOP) 
or formulation (4) (for differentiable MOP), and solve it. 


Step 5- Present the obtained (weakly, properly) efficient solution(s) and the corresponding 


objective function vector(s) to the DM. Let DM chooses one of them. In this case, 
different states can occur: 


5.1. Ifthe DM approves this solution as the most preferred one, denote this solution 
by x*+1 and go to Step 6. 


5.2. If the DM wants to obtain other solutions with the same preference informa- 
tion, go to Step 4. Note that, in this case, Step 4 should be executed with 
other values for parameters and index sets. 


5.3. If the DM wants to provide new preference information, denote this solution 
by x*+1, set k :=k+1 and go to Step 3. 


Step 6- Stop. 
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corresponding scalarizing formulation (3) has a reduced number of con- 
straints which causes a decrease in solving time. 


Unlike the algorithms proposed in [21,28,32], the GIA allows the DM to 
specify his/her preference information in six different ways. Since the 
satisfaction of DM is an important factor in the interactive algorithms, 
this aspect of the GIA will increase the satisfaction of DM. 


To propose a general algorithm for solving an MOP, it is necessary to 
convert the MOP problem to a general scalarized problem with per- 
haps some additional constraints. It is obvious that the number of 
constraints added to the general scalarized problem has a major effect 
on the computational time. In Table 1 (for nondifferentiable MOPs) 
we compare the number of constraints added to the suggested gen- 
eral formulation (3) with those added to some general formulations as 
GLIDE [21], GLIDE-II [28] and GENWS [32]. 


Table 1: Number of additional constraints in each formulation (in nondiffer- 
entiable case) 


Methods GENWS GLIDE GLIDE-II nondiff Our formulation 


Reference direction approach - - 0 0 


GUESS card(J5 ) 2Qp 0 0 


STOM P 
SPOT - 
Modified reference point card( Js) + card(J®) 
RD - 
€— constraint 
Weighted sum 
Hybrid 
ISWT 
PROJECT 
Weighted Tchebycheff 


0 0 
p-1 
card(J®) card(J®) 
- card(J$) + card(J$) 
= p= 1 
0 0 


iS) 
as] 
3 
| 
BH 


brit 
ly 


(d) 


One of the most important theoretical advantages of the proposed gen- 
eral formulations is that, Theorem 5 enables us to provide results con- 
cerning proper efficiency. Unboundedness of the trade-offs means, prac- 
tically, ignoring at least one of the objective functions when the DM 
wants to improve another objective function, which is not satisfactory 
to the DM. Since properly efficient solutions have bounded trade-offs, 
the DM can improve some unacceptable objective functions with no 
concern. 


All the provided theorems where established without convexity assump- 
tions. In fact, the main MOP can be convex or nonconvex. 


7 A numerical example 


In this section, we illustrate the procedure mentioned in the GIA on an 
engineering example of designing a four-bar plane truss, studied in [5]. This 
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problem, has two conflicting objective functions. We should minimize the 
volume of the truss (f1), and its joint displacement (f2), subject to given 
physical restrictions on the feasible cross-sectional areas 71, 22,23, and x4 of 
the four bars. The stress on the truss structure is caused by several forces of 
magnitude F, and 2F. The length L of each bar and the elasticity constants EF 
and o of the materials involved are modelled as constants. The mathematical 
model of this example is as follows: 


Minimize fi(x) = L(2a, + Varo + V2r3 + 24), 


FL,2 #2 2/2 1 
fhe) = a(S + v ve aay 
F< <¥, 
ot UV AG) Sm SF, 
~ | V2(9)<a3 < %, 
Peay < ¥, 
where, the constant parameters are chosen as F = 10 kN, E = 2 x 


10° kN/cm?, L = 200 cm and o = 10 kN/cm?. 


The ideal and nadir values for objective functions of this problem are 
obtained as y! = (y/,y4) = (1400, -—5.7191 x 1074) and yN = (y¥,y’) = 
(3.4971 x 103,0.0406). Now, based on the GIA, at first, we should find an 
initial solution. To this end, the «—constraint scalarizing problem (13) is 
used, which can be obtained from the proposed formulations by parameters 
and index sets given in Subsection 4.7, with | = 2 and €; = 1800. By solving 
the obtained problem, we find (1.3906, 1.9963, 1.4142, 1.3957) for variables, 
and (1800,0.0157) for objective functions. As it can be seen, the values 
of the objective functions are between ideal and nadir values. Let #1 = 
(1.3906, 1.9963, 1.4142, 1.3957) and f(#!) = (1800, 0.0157) are shown to DM. 


Suppose, the DM wishes to express his/her preference information as the 
reference point y' = (yt, y3) = (1600, 0.01). Based on this preference given 
by DM, one of the reference point based scalarizing problems can be selected. 
Here, we set parameters in our formulation, such that the GUESS scalariz- 
ing problem is obtained, and by solving it, (1.4613, 2.0666, 1.4142, 1.4613) is 
obtained for variables and (1861.3, 0.0142) is attained for the objective func- 
tion values. At this iteration, the volume of the truss has increased and its 
joint displacement has slightly decreased. Now, assume that the DM wants 
to change the type of his/her preference information. According to Step 5 of 
the GIA, set #7 = (1.4613, 2.0666, 1.4142, 1.4613), f(#7) = (1861.3, 0.0142) 
and k = 2. Now, Step 3 is executed. 


Assume that the DM wants to classify the objective functions in two 
classes J? = {1}, and J? = {2}. This means, the DM wants to improve f; by 
somewhat relaxing f2. Assume that the DM gives us ¥? = (1500, 0.03). Now, 
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the RD scalarizing problem is used with a = 0.5. By solving the problem, 
(1.1876, 1.6796, 1.4142, 1.1876) is obtained for variables and (1587.6, 0.0221) 
is attained for the objective functions. Assume that the DM wants to provide 
new preferences by selecting weights for the objective functions. To this 
end, let 2° = (1.1876, 1.6796, 1.4142, 1.1876), (#3) = (1587.6, 0.0221), and 
k = 3. Then, Step 4 is executed by (w}, w3) = (3,4) as weights given by 
the DM. By solving the weighted sum problem (14), (1, 1.4142, 1.4142, 1) is 
obtained for variables and (1400, 0.03) is obtained for the objective values. As 
it can be seen, the volume of the truss is in its ideal value, and this satisfies the 
DM. It is important to point out that by Theorem 5, the obtained objective 
vector is a properly nondominated point. 


8 Conclusions 


In this article, we suggested a general scalarizing formulation to obtain a 
global interactive algorithm for multiobjective optimization problems. We 
proposed the formulation in two versions; one of them for differentiable and 
the other for nondifferentiable MOPs. By selecting suitable values for pa- 
rameters, we proved that optimal solutions of the suggested general scalar- 
izing problem are (weakly, properly) efficient solutions for the main mul- 
tiobjective problem. Moreover, it was shown that many scalarizing prob- 
lems used in different interactive methods as GUESS, reference direction 
approach, Step, STOM, SPOT, modified reference point, and RD methods 
can be obtained from the proposed general formulation, by selecting suitable 
transformations. Some of the popular scalarizing problems such as, weighted 
sum, €—constraint, and hybrid problems derived from our general scalarizing 
problem. In addition, we proposed a general interactive algorithm. In the 
proposed algorithm, the DM could express his/her preference information in 
six different ways, and based on the kind of information given by the DM, 
a suitable scalarizing problem, by selecting appropriate values for parame- 
ters and index sets in the general formulation, was selected. Finally, by a 
numerical example we illustrated that how the proposed general interactive 
algorithm can be used. 


However, for the future investigation, developing a software based on the 
suggested general interactive algorithm can be worthwhile. Also, proposing a 
general interactive procedure for approximate efficient solutions of an MOP 
can be worth studying. To this end, studying three recently published papers 
by Ghaznavi-ghosoni! and Khorram [10], Ghaznavi-ghosoni et al. [11] and 
Ghaznavi [12] is recommended. 


1 Previous name of the first author 
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